Transverse Commensurability Effect For Vortices in Periodic Pinning Arrays 
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Using computer simulations, we demonstrate a new type of commensurability that occurs for 
vortices moving longitudinally through periodic pinning arrays in the presence of an additional 
transverse driving force. As a function of vortex density, there is a series of broad maxima in 
the transverse critical depinning force that do not fall at the matching fields where the number of 
vortices equals an integer multiple of the number of pinning sites. The commensurability effects are 
associated with dynamical states in which evenly spaced structures consisting of one or more moving 
rows of vortices form between rows of pinning sites. Remarkably, the critical transverse depinning 
force can be more than an order of magnitude larger than the longitudinal depinning force. 

PACS numbers: 74.25.Qt 



C ■ 

o ■ 
o , 
i 

s : 

a ■ 



c 

o 
o 



> 

ON 

in 

OO 

o 

OO 

o 



X 



Matching effects for vortices in periodic pinning ar- 
rays have been studied extensively for different types of 
pinning lattice geometries [H, 0, 0, 0, 0, 0] • As a func- 
tion of magnetic field, the critical current passes through 
a series of peaks generated by commensurability effects 
that occur when the number of vortices equals an integer 
multiple of the number of pinning sites and the vortex 
ground state is an ordered crystalline structure [1, 0, 0] . 
When each pinning site can capture only one vortex, the 
excess vortices at fields above the first matching field are 
located in interstitial sites and depin first under an ex- 
ternal drive 0, 0, 0, 0, 0] • Once the interstitial vortices 
are moving under a longitudinal drive, it is possible to 
apply an additional transverse drive in the direction per- 
pendicular to the vortex motion. In this case, although 
the vortices are mobile in the longitudinal direction, they 
can remain pinned in the transverse direction and there 
can be a finite transverse critical depinning threshold. 

The possibility of a transverse depinning threshold for 
moving vortices was initially predicted for systems with 
random pinning when the moving vortices form well de- 
fined channels [8( , and transverse depinning thresholds in 
randomly pinned systems have been observed in numeri- 
cal simulations Q and experiments 1(J 11 1 . For vortices 
moving in the presence of a periodic pinning array, a fi- 
nite transverse depinning threshold has been measured 
at high drives when all of the vortices are moving and 
has a value that depends on the angle between the lon- 
gitudinal driving direction and a symmetry axis of the 
pinning lattice. Here, the most prominent transverse de- 
pinning thresholds and dynamical locking effects occur 
for driving along the principal axes of the pinning lat- 
tice [13, El ■ This type of effect has been experimentally 
observed for colloidal particles moving over periodic sub- 



strates 14J|. When only one vortex can be captured by 



each pinning site, motion of the vortices at low drives oc- 
curs as a flow of interstitial vortices between vortices that 
remain trapped at the pinning sites [2|, [3|, [6( ■ In this case, 
it is not known whether a transverse depinning threshold 
exists, and in general it is not known how the transverse 



depinning threshold varies with magnetic field. 

It might be expected that the transverse depinning 
threshold would simply exhibit peaks at the same mag- 
netic fields where peaks in the longitudinal depinning 
threshold appear. In this work, we demonstrate that al- 
though there are enhancements of the transverse depin- 
ning threshold at certain fields, these fields are not related 
to fields which produce peaks in the longitudinal depin- 
ning threshold, but are instead associated with dynamical 
matching conditions. The distinct dynamical matching 
effects appear because the moving vortices assume a dif- 
ferent structure than the static vortex ground state. Dy- 
namical commensurability effects occur when an integer 
number of moving interstitial vortex rows form between 
adjacent rows of pinning sites. The dynamical matching 
effects are much broader than the static matching effects 
and have maxima that encompass several static match- 
ing fields. An oscillatory critical current appears for the 
dynamical transverse commensurability effect. This is 
similar to the critical current oscillations seen for vor- 
tices depinning in artificial channels (lEL flit ] or in lay- 
ered or strip geometries [121,1131, although in the channel, 
layer, or strip systems, the commensurations arise due to 
matching effects of the vortex ground state rather than 
the dynamical matching effects observed in the present 
work. Remarkably, we find that the transverse depinning 
threshold can be up to an order of magnitude larger than 
the longitudinal depinning threshold. 

We numerically simulate a two-dimensional system 
with periodic boundary conditions in the x and y di- 
rections containing N v vortices and N p pinning sites fol- 
lowing a procedure similar to that used in previous sim- 
ulations for vortices in periodic pinning arrays. The 
number of vortices is proportional to the applied mag- 
netic field B = Bz, which is normal to our simulation 
plane. The repulsive vortex-vortex interaction force is 
given by F™ = E$" ? - /o^i (#y A)Ry , where K x is a 
modified Bessel function, Bij = |Rj — Rj| is the dis- 
tance between vortex i and j located at R; and Rj, 
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(Ri-RfV-Rtf, f = ^/(2^/i A 3 ), = h/2e is the 



2 



elementary flux quantum, and A is the London penetra- 
tion depth. The pinning sites are placed in a triangular 
lattice, and the field at which the number of pinning sites 
equals the number of vortices, N p = N v , is defined as the 
matching field B^. The individual pinning sites are mod- 
eled as parabolic traps of radius r p — 0.35A and strength 
FP = 1.25, withFf = £f» FPf (R lk /r p )O(r p ~R lk )R tk , 
where 8 is the Heaviside step function, Ri k = |Rj — R?|, 
Rife = (Ri — R%)/Rij, and R| is the location of pin k. 
The overdamped equation of motion for a single vortex i 
is 

ait, = p p Feat 

where 77 = 1 is the viscous damping term. F ea:t repre- 
sents the net force from an applied current and is given 
by F ext = F^f x + F^ r f y, where the longitudinal drive 
is applied in the x direction and the transverse drive 
F^ r is applied in the y direction. The initial vortex po- 
sitions are obtained by simulated annealing. The drive is 
first applied in the longitudinal direction in increments of 
AF^ = 0.0015, with 15000 simulation time steps spent 
at each current increment. Once the longitudinal drive 
reaches the desired value, it is held fixed while the trans- 
verse drive is increased from zero with the same current 
increment protocol. We find that our increment rate is 
sufficiently slow to avoid any transient effects. The longi- 
tudinal and transverse critical depinning thresholds, 
and F^ T , are obtained by measuring the vortex velocity 
(V a ) — N^ 1 Yli" v » ' oti with a = x,y, and identifying 
the driving force at which (V a ) > 0.001. 

We first study a system with a low pinning density to 
ensure that a portion of the vortices are located in the 
interstitial sites. The existence of a clearly defined depin- 
ning threshold which varies nonmonotonically with field 
is illustrated in Fig. 1(c), where we show {Vr r ) versus 
Fl r for B/B^ = 4.33, 6.0, 8.0, and 10.0 in a system 
with F = O.O520 O /A 2 and fixed Fp = 0.6. From a se- 
ries of simulations, we obtain the variation in F^ r versus 
BjB^ plotted in Fig. 1(a). Four well defined maxima 
in F^ r appear that are centered near B/B^ = 2.0, 6.0, 
12.0, and 17.0. Figure 1(b) shows a blowup of the region 
B/B^ > 2.0, where the oscillation in F ( 7 V can be seen 
more clearly. This oscillation is distinct from the match- 
ing effects observed for longitudinal depinning P, H, 0, , 
where well defined peaks occur at integer matching fields. 
The maxima in Fig. 1(a) are much broader than in the 
longitudinal depinning case and each encompass three or 
more matching fields. Similarly, the minima in F^ r also 
each spread over several values of B/B^. 

For all fields BjB^ > 1.0, we find that F c Tr is signif- 
icantly larger than the longitudinal critical force Fjf, as 
shown in Fig. 1(d) where we plot (Vl)(N v /N p ) versus Fp 
and (Vrr)(N v /N p ) vs F^ r . Here the velocities have been 
scaled by N p rather than N v for presentation purposes. 
For B^ = 2.0, the transverse depinning threshold F^ r is 




FIG. 1: (a) The transverse critical depinning force F^ r vs 
B/B^ for a system with B<j> = 0.052<j!>o/A 2 and fixed longi- 
tudinal drive = 0.6. The points a, c, and e refer to the 
fields illustrated in Fig. 2. The maxima are labeled according 
to the number of moving vortex rows between adjacent pin- 
ning rows: 1R, one row; 2R, two rows; 3R, three rows; and 
4R, four rows, (b) A blow up of panel (a) for BjB^, > 2.0. 
The points a, c, and e refer to the fields illustrated in Fig. 3. 
(c) The transverse velocity (Vrr) vs transverse force F^ r for 
B/B$ = 4.33 (open triangles), 6.0 (filled circles), 8.0 (filled 
squares), and 10 (open diamonds), (d) The scaled longi- 
tudinal velocity (Vl)(N v /N p ) (filled circles) versus longitu- 
dinal drive F^ and transverse velocity (VTr)(N v /N p ) (open 
squares) vs transverse drive F£ r for B /B^ — 6.0 (left curves) 
and 2.0 (right curves). 



about six times higher than the longitudinal depinning 
threshold F C L . Both depinning thresholds are lower for 
B/B^ = 6.0; however, Fj r is again much higher than 
F C L - 

In Fig. 2(a) we show the vortex and pinning site po- 
sitions for point a in Fig. 1(a) at B/B^ = 2.67 and 
in Fig. 2(b) we illustrate the vortex trajectories for 
Fp r < F^ r , just below the transverse depinning transi- 
tion. There is a single row of moving interstitial vortices 
between neighboring rows of pinning sites and the vortex 
lattice is anisotropic, with higher vortex density in the in- 
terstitial rows than in the pinned rows. The same vortex 
structure appears for 1.0 < B/B^ < 2.9, correspond- 
ing to the maximum in F^ r marked 1R in Fig. 1(a). In 
Fig. 2(c), we plot the vortex positions for B/B^ = 4.0 
at a minimum of F^ r found at the point marked c in 
Fig. 1(a). The interstitial rows are no longer uniform and 
consist of an interlacing of double rows with single rows. 
Figure 2(d) shows that the vortex trajectories at this field 
are more disordered. The vortex positions and trajecto- 
ries at point e in the region marked 2R in Fig. 1(a) for 
B/B^ = 6.0 appear in Fig. 2(e,f). Here there are two 
well defined rows of moving vortices between adjacent 
pinning site rows. 
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FIG. 2: The vortex positions (black dots), pinning site loca- 
tions (open circles), and vortex trajectories (black lines) for 
the system in Fig. 1(a). (a), (b) BjB$ — 2.67, marked a in 
Fig. 1(a). (c), (d) B/Bj, = 4.0, marked c in Fig. 1(a). (e), (f) 
B/B<fy = 6.0, marked e in Fig. 1(a). 




FIG. 3: The vortex positions (black dots), pinning site loca- 
tions (open circles), and vortex trajectories (black lines) for 
the system in Fig. 1(a). (a), (b) B/B^ = 8.0, marked a in 
Fig. 1(b). (c), (d) B/Bj, = 10.67, marked c in Fig. 1(b). (e), 
(f) BjB^ = 17, marked e in Fig. 1(b). 

We find that maxima in F^ r occur whenever there 
is an integer number of moving rows of interstitial vor- 
tices between neighboring pinning rows. Since the num- 
ber of vortices in each interstitial row can vary over a 
considerable range without destroying the row structure, 
the maxima in F^F r are much broader than the peaks in 
associated with longitudinal commensuration effects. 
The row structures become increasingly anisotropic with 
increasing field until a buckling transition occurs which 
marks the end of the maximum in F<F r . In Fig. 3(a) we 
illustrate the vortex positions for B/B^ = 8.0 at a min- 
imum in F^F r found at the point marked a in Fig. 1(b). 
The interstitial vortices form a mixture of two and three 
interstitial rows between pinning site rows, producing the 




FIG. 4: (a) The longitudinal critical depinning force F C L vs 
for B/Bff, — 2.0 (filled squares) and B/B^ = 2.5 (open 
squares) and the transverse critical depinning force F^ r vs 
B<f> for B/B$ = 2.0 (filled circles) and B/B^ = 2.5 (open 
circles) (b) F^ r vs the applied longitudinal force F^/F p for 
B/B4, = 3.0 (filled circles) and B/B^ = 3.67 (filled squares). 

nonuniform trajectories shown in Fig. 3(b). At the max- 
imum in F^F r marked c in Fig. 1(b), corresponding to 
B/B^ = 10.67, Fig. 3(c,d) shows that there are three well 
defined rows of moving vortices between adjacent pinning 
site rows. Similarly, Fig. 3(e,f) indicates that there are 
four interstitial vortex rows at B/B^ = 17, which falls on 
the maximum in F^ r at the point marked e in Fig. 1(b). 
Near B/B$ = 14, where Fj r passes through a minimum, 
the interstitial vortices form a mixture of three and four 
rows, while for B/B^ > 19 there is a mixture of four and 
five interstitial rows (not shown). 

Commensurability effects generated by the presence of 
an integer number of vortex rows between line-like bar- 
riers have been observed for longitudinal vortex motion 
through channel geometries flBL as well as critical 
currents in layered materials [171 ] , superconducting strips 
[lij ]. and anisotropic pinning arrays [jjj . In all these 
cases the commensurability occurs in the static vortex 
configurations. This is distinct from the transverse de- 
pinning maxima that we observe here, which arises due 
to commensurations in the dynamical interstitial vortex 
configuration. 

As shown in Fig. 1(c), the transverse depinning thresh- 
old is much higher than the longitudinal depinning 
threshold. In Fig. 4(a) we quantify this effect by plot- 
ting F^ and F^F r as a function of pinning density B^ 
for a commensurate field B/B$ = 2.0 where F^ passes 
through a peak and for an incommensurate field B / B^ = 
2.5. At the incommensurate field, F C L and F^ r are both 
reduced. At B/B<f, = 2.0, F C L increases monotonically 
with B^ while F^ r shows a smaller increase; however, 
F^ r is significantly larger than F^ over the entire range 
of pinning densities studied. At the incommensurate field 
B/Bff, = 2.5, we find a similar trend; however, F C L in- 
creases much more slowly than FjF r with increasing B^ 
and at B^ = 0.6, F^F r is nearly an order of magnitude 
larger than F C L . In addition to increasing with increasing 
B^, the ratio Fj r '/F£ increases with decreasing B^ as 
B^ approaches zero due to the different rates at which 
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the two thresholds approach zero. 

The higher value of F^ r compared to F C L can be under- 
stood by considering that longitudinal dcpinning occurs 
from the ground state configurations of the interstitial 
vortices 0, 01 and is determined by the repulsive inter- 
actions between the vortices at the pinning sites and the 
interstitial vortices. In the ground state, the interstitial 
vortices occupy positions that lower the repulsion from 
the pinned vortices, and the initial longitudinal depin- 
ning occurs when the interstitial vortices begin to move 
between the pinning sites, such as in Fig. 2(a). For the 
transverse depinning, when the interstitial vortices are 
moving at a sufficiently high velocity in the longitudi- 
nal direction they do not have time to slip between the 
pinned vortices in the transverse direction, but instead 
come into close proximity with the pinned vortices and 
interact strongly with them, resulting in a high repul- 
sive barrier for depinning. If the longitudinal drive is set 
to a lower value before the transverse drive is applied, 
the interstitial vortices have more time to pass between 
the pinned vortices and F^ r decreases. In Fig. 4(b) we 
illustrate this effect by plotting Fj r versus Fj^/F p for 
B/B^ = 3.0 and B/B^ = 3.67. In both cases F c Tr 
increases from a low value with increasing F^/F p un- 
til reaching a maximum value at Fg/F p = 1.07 for 
B/Bcj, = 3.0 and at F%/F p = 1.12 for B/B^ = 3.67. 
Above this drive, Fj r decreases with increasing F^/F p 
as Ff) I F p approaches 1 since the vortices at the pinning 
sites begin to depin for the higher longitudinal drives, re- 
ducing the magnitude of the transverse critical current. 

In summary, we have shown that a new type of dy- 
namical commensurability effect can occur for vortices 
in periodic pinning arrays. When interstitial vortices are 
moving between pinned vortices and an additional trans- 
verse force is applied, there is a finite transverse critical 
depinning force which oscillates with field. The oscilla- 
tion is not simply related to the matching of the vortices 
with the number of pinning sites as in the case for the lon- 
gitudinal depinning, but is associated with the dynamical 
structure of the vortices which allows for integer or non- 
integer numbers of rows of moving interstitial vortices 
between adjacent rows of pinning sites. The transverse 
commensurability effects are much broader than those 
seen for the longitudinal depinning and each maximum 
in the transverse depinning force can span several match- 
ing fields. Remarkably, the transverse depinning force 
can be more than an order of magnitude larger than the 
longitudinal depinning force due to the fact that the mov- 
ing interstitial vortices are unable to move between the 
pinned vortices without coming close to the pinned vor- 
tices, which creates a strong repulsive barrier for trans- 
verse depinning. 

This work was carried out under the auspices of the 
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